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Abstract 

MicroRNAs(miRNAs) are transcriptional and posttranscriptional regulators involved in nearly all known biological processes in distant 
eukaryotic clades. Their discovery and functional characterization have broadened our understanding of biological regulatory mech- 
anisms in animals and plants. They show both evolutionary conserved and unique features across Metazoa. Here, we present the 
current status of the knowledge about the role of miRNAin development, growth, and physiology of teleost fishes, in comparison to 
other vertebrates. Infraclass Teleostei is the most abundant group among vertebrate lineage. Fish are an important component of 
aquatic ecosystems and human life, being the prolific source of animal proteins worldwide and a vertebrate model for biomedical 
research. We review miRNA biogenesis, regulation, modifications, and mechanisms of action. Specific sections are devoted to the role 
of miRNA in teleost development, organogenesis, tissue differentiation, growth, regeneration, reproduction, endocrine system, and 
responses to environmental stimuli. Each section discusses gaps in the current knowledge and pinpoints the future directions of 
research on miRNA in teleosts. 
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Introduction 

Small nonprotein-coding RNAs (ncRNAs) are short, 18-40 nu- 
cleotide (nt) sequences, with diverse biogenesis pathways and 
regulatory mechanisms. They have convoluted relationships, 
in which they cooperate, compete, or regulate each other 
(Ghildiyal and Zamore 2009). They are involved in basic cellular 
processes, including differentiation, proliferation, and apopto- 
sis (Bartel 2009). 

Among several types of small RNAs, microRNAs (miRNAs), 
approximately 22-nt-long posttranscriptional regulators of 
mRNA, have been intensively investigated in recent years. 
They were discovered in early 1990s (Lee et al. 1993; 
Wightman et al. 1993), but it took almost a decade to discover 
how massive their involvement in gene expression regulatory 
networks is. It is estimated that in mammals over 60% of 
mRNAs have conserved miRNA-binding sites (Friedman et al. 
2009; Guo et al. 2010). Also in teleost fishes, miRNAs are 
involved in the development and various physiological pro- 
cesses (Wienholds et al. 2005; Salem et al. 2010; Bizuayehu 
et al. 2012a; Mishima 2012; Wei et al. 2012). Although no 
estimation or quantification has been made yet to determine 
the extent of miRNA participation in regulatory network in 
teleosts, it can be anticipated that a considerable set of 
mRNAs is under their modulation, based on high conserva- 
tion of miRNA among animals in general and vertebrates 



in particular. Teleost miRNAs were first reported in zebrafish 
(Lim et al. 2003) and miRNA repositories during zebrafish de- 
velopment, and some functions were characterized (Chen 
et al. 2005; Giraldez et al. 2005; Wienholds et al. 2003, 
2005). Although a number of studies on miRNA in some 
other teleosts has been reported since then (Salem et al. 
2010; Fu et al. 2011; Xia et al. 2011; Barozai 2012; 
Andreassen et al. 2013; Bekaert et al. 2013; Xu et al. 2013), 
information on the functions of miRNAs in teleosts has largely 
been obtained from studies carried out on zebrafish. 

Teleosts are the most speciose among the vertebrate line- 
age with an estimated number of species exceeding 25,000 
(Nelson 2006); therefore, miRNA characterization and func- 
tional studies performed in few species so far are presumptu- 
ous to conclude on general principles. Currently, there are 
1,250 miRNAs identified in 8 teleost species, representing 5 
orders (miRBase v. 20), which reflects how little has been done 
in characterization of miRNAs in fishes when compared with 
mammals (table 1). 

Fish are important in broad terms of ecology and food 
production. Both exploitation and conservation tasks need a 
baseline knowledge of the habitat and physiology of a spe- 
cies in question. Understanding molecular mechanisms and 
functions can provide sustainable and more efficient, knowl- 
edge-based solutions. Teleosts also serve as biological models. 
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Genetic studies in teleosts have widened our knowledge of 
number of biological pathways that are common in verte- 
brates (Howe et al. 2013). Teleosts show a number of advan- 
tages, including high fecundity, oviparity, quick development, 
easiness of manipulation and production of genetic modifica- 
tions, early development of functional systems, and the po- 
tential of tissue regeneration. Studies of teleost genetic 



Table 1 

Number of Precursor and Mature miRNAs Identified in Teleost Fish 
and Compared with Other Vertebrates (miRBase v.20) 





Precursor 


Mature 


Fish 


1,250 


1,044 


Cyprinus carpio 


134 


146 


Danio rerio 


346 


255 


Fugu rubripes 


129 


108 


Hippoglossus hippoglossus 


40 


37 


Ictalurus punctatus 


281 


205 


Oryzias latipes 


168 


146 


Paralichthys olivaceus 


20 


38 


Tetraodon nigroviridis 


132 


109 


Amphibians 


211 


196 


Reptiles 


282 


416 


Birds 


980 


1,330 


Mammals 


9,076 


11,717 



regulatory elements, such as miRNAs, can provide a much 
needed insight into the human gene regulatory networks 
through orthologous gene functional studies. Given the fact 
that 82% of human genes that are implicated in genetic-re- 
lated diseases have their orthologs in zebrafish (Howe et al. 
201 3), many challenges in human medicine can be addressed 
by better understanding of conserved genes and molecular 
mechanisms. 

The aim of this review is to summarize the recent progress 
made in teleost miRNA research and discuss areas of future 
studies. We review the general features of miRNAs, their roles 
in teleost development and physiology, and give an overview 
of modeling miRNA functions in teleosts. We discuss the 
major gaps in knowledge on miRNA in teleosts, particularly 
in regard to other model systems. 

Overview of miRNA Biology 

Biogenesis and Mechanisms of Action 

There are several known pathways of miRNA biogenesis 
(fig. 1). In the canonical pathway, miRNA synthesis begins in 
the nucleus, where miRNA genes are transcribed by RNA 
polymerase II (or polymerase III for some miRNAs) and form 
capped and polyadenylated primary nascent transcripts 
(pri-miRNAs) of variable length, ranging from hundreds to 



Cytoplasm 



a ' Canonical pathway b ' Alternative pathway 




Fig. 1. — miRNA biogenesis pathways. miRNA processing starts from the nucleus, (a) Canonical pathway represents those miRNAs that are transcribed by 
polymerase II and then are processed by Drosha and associated proteins, (b) Alternative miRNA processing pathway represents those miRNAs that bypass 
Drosha processing; for example, lariats that debranch and form pre-miRNA structure, (c) Pre-miRNA processing in cytoplasm by Dicer together with other 
RNA-binding proteins, (d) miRISC binds to 3'-UTR of an mRNA for translational repression, (e) miRISC is transported to the nucleus and binds to ncRNAs 
including pri-miRNA to repress their processing or interfere with their functions. 
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thousands ribonucleotides (Cai et al. 2004; Borchert et al. 
2006). These transcripts can be monocistronic (single hairpin) 
or polycistronic (multiple hairpins). A hairpin has three defined 
features: Terminal loop, internal bulges, and double-stranded 
stem. 

The pri-miRNAs are further processed to shorter (-70 nt) 
precursor miRNAs (pre-miRNAs) by Drosha (Lee et al. 2003), 
an RNase type III enzyme, together with at least 20 other 
polypeptides, such as DiGeorge syndrome critical region 8 
(Dgcr8, also termed as Pasha) and a double-stranded RNA- 
binding domain protein (Tomari and Zamore 2005) (fig. 1a). 
The resulting hairpin structure has phosphate and hydroxyl 
groups at its 5'- and 3'-ends, respectively, and a characteristic 
2-nucleotide overhang at the 3'-end (Han et al. 2004; Tomari 
and Zamore 2005; Kim et al. 2009). Pre-miRNAs are regulated 
by diverse processes (Burroughs et al. 2011). Some miRNAs 
require additional protein factors, such as p68 or p53 (Fukuda 
et al. 2007) and coactivator KH splicing regulatory protein 
(Trabucchi et al. 2009). However, pre-miRNAs can also be 
formed through Drosha-independent pathway (fig. 1b) 
(Okamura et al. 2007; Ruby et al. 2007). It is not clear, 
when Drosha-dependent mechanism emerged in metazoan 
evolution. Nevertheless, mirtrons, that is miRNAs that are pro- 
cessed from mRNA introns by a spliceosome, do exist in dis- 
tant animal lineages, such as nematodes (Ruby et al. 2007), 
insects (Okamura et al. 2007), and mammals (Babiarz et al. 
2008), as well as in plants (Zhu et al. 2008). 

After the transportation of a pre-miRNA from nucleus to 
cytoplasm, which is facilitated by Exportin-5 in the presence of 
Ran-GTP (Yi et al. 2003; Bartel 2004; Zeng and Cullen 2004; 
Kim et al. 2009), pre-miRNA is diced by Dicer, an another 
RNase III, and associated proteins, such as trans-activator 
RNA-binding protein, protein activator of PKR (PACT), and 
Argonaute 2 (Ago2); this process yields in approximately 
22-nt-long miRNA duplex (Lee et al. 2002, 2003). However, 
some miRNAs origin from a Dicer-independent pathway 
(Chelouf i et al. 201 0; Dueck and Meister 201 0). In the canon- 
ical miRNA processing pathway, one of the strands is loaded 
to Ago protein through a mechanism unresolved yet (fig. 1c). 
Several hypothetical models explain the incorporation of a 
single strand to form an active miRNA-induced silencing 
complex (miRISC), including active incorporation of one of 
the strands after unwinding by a helicase using ATP 
(Salzman et al. 2007), ATP-dependent loading of miRNA 
duplex to Ago but passive process of strand selection 
(Kawamata et al. 2009) and stepwise loading, wedging, and 
unwinding by Ago (Kwak and Tomari 2012). 

The major role of miRNAs in cellular processes is posttran- 
scriptional repression of mRNA in cytoplasm. However, recent 
studies indicate that mature miRNA can be imported into the 
nucleus and repress ncRNAs (Leucci et al. 2013). The posttran- 
scriptional suppression is predominantly achieved by binding 
the miRISC at 3'-UTR of an mRNA. Base complementarity be- 
tween miRNA and mRNA influences the final outcome of the 



repression (fig. 1d), in which a perfect base pairing results in 
target degradation, whereas imperfect base pairing yields se- 
questration of a target (Bartel 2009). The "seed" sequence, 
nucleotides at the positions 2-8 from the 5'-end of a mature 
miRNA, is the major determinant of imperfect matching. Seed 
sequences are highly conserved among species and used to 
categorize different miRNA families (Griffiths-Jones et al. 
2008). Different models have been proposed to explain the 
mechanism of interaction between miRISC and mRNA, includ- 
ing seed matching, which is the complementarity between the 
seed and its target mRNA (Lee et al. 1993), seedless matching 
(Lai et al. 2009), centered paired site (Shin et al. 2010), and 
pivot pairing and transitional nucleation models (Chi et al. 
2012) (fig. 2). In mammals, this interaction can result in 
reduced translational initiation rate followed by mRNA dead- 
enylation and degradation (Guo et al. 2010; Hu and Coller 
2012). Different mechanisms of miRNA-mediated gene 
silencing have been described (Fabian and Sonenberg 2012). 

Although the most known function is to mitigate mRNA 
translation, animal miRNAs are involved also intranslational 
promotion and repression of other noncoding RNAs. 
Posttranscriptional repression by miRNAs can be achieved 
through inhibition of translation initiation, inhibition of trans- 
lation elongation, premature termination of translation, or 
deadenylation (Eulalio et al. 2008; Moretti et al. 2012). 
In zebrafish, translational control by miRISC is largely observed 
before gastrulation; however, during and after the gastrula- 
tion, the action of miRISC is switched to mRNA destabilization 
(Subtelny et al. 2014). miRISC binds poly(A)-binding protein 
and recruits deadenylase to promote poly(A)-tail cutting 
(Fabian et al. 2009) (fig. 3a). Pasquinelli and Ruvkun (2002) 
suggested that miRNAs can bind to the 5'-end of mRNAs, pre- 
mRNAs, and DNA to suppress translation, facilitate alternative 
splicing, and form RNA-DNA duplexes, respectively. miRNAs 
bind not only to untranslated regions (UTRs) but also to exons 
(fig. 3b) (Tay et al. 2008; Fang and Rajewsky 2011). 
Furthermore, miRNAs can be imported to nucleus and repress 



miRNA mRNA ,S> 
Perfect match 5' nm A(n)3' m> / 
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-A(n)3' 
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Fig. 2. — Four miRNA-mRNA interaction models and their final out- 
come. A perfect pairing between miRNA (green) and mRNA (blue) results 
in the degradation of mRNA, whereas imperfect matching (only seed, 
seedless, and transitional nucleation pairings) results in translational 
inhibition. 
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3' UTR 
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Fig. 3. — Examples of miRNA mechanisms of action, (a) mRNA translation includes initiation, elongation, and termination steps, which are facilitated by 
the binding of different RNA-binding proteins. Poly(A)-binding protein (PABP) binds to poly(A) tract that in turn binds to elF4G. elF4G serves as a platform for 
the binding of elF4E (binds to m 7 G cap structure at the 5'-end of mRNA), elF4A, elF3, and other proteins. These interactions shape mRNA and enhance the 
translation (left). However, the interaction of PABP with miRISC augments miRNA-mediated translational repression through the recruitment of deadenylase 
(right, top). In absence of PABP, miRISC binds to 3'-UTR and destabilizes an mRNA (right, bottom), (b) miRISC binds to exons for translational repression 
possibly by limiting translational elongation. The scheme is simplified and depicts only some of the RNA-binding proteins taking part in the process. 



other ncRNAs (fig. 1e). For instance, MALAT1, a long ncRNA, 
is a target of miR-9 in the nucleus (Leucci et al. 2013). 
Moreover, miRNAs have functions in translational promotion; 
for example, miR-1 0 is implied in enhancing mRNA translation 
of a ribosomal protein by binding at the 5'-UTR (0rom et al. 
2008). These different mechanisms highlight the flexibility of 
miRNA action. 

Regulation 

Strict c/s- and frans-acting regulatory mechanisms exist in a cell 
to control miRNA biogenesis at different levels. These regula- 
tory steps can be categorized as transcriptional regulation, 
posttranscriptional nuclear regulation (microprocessor, shut- 
tle, or autoregulation), posttranscriptional cytosolic regulation 
(processors regulation, loading regulation, or strand selection), 
and decay (fig. 4). 

Transcriptional Regulation 

miRNA genomic location can be intergenic, intronic, inside 
repetitive elements, or a standalone gene with its own pro- 
moter (fig. 4a). Genomic location of a miRNA determines its 
transcription, which depends on promoter and enhancer ele- 
ments (Cai et al. 2004; Borchert et al. 2006). A number of 
miRNAs is clustered. Some miRNA clusters have multiple 
miRNA promoters inside a cluster, such as C19MC in mam- 
mals (Bortolin-Cavaille et al. 2009). Some miRNAs are regu- 
lated together with their targets, such as miR-1 0c and HoxB4a 



in zebrafish (Woltering and Durston 2008). A miRNA and its 
target can be transcribed as a single transcriptional unit, for 
example, miR-26b and ctdsp2 (Han et al. 2012), or miR-412 
and Mirg (Melamed et al. 2013). 

Posttranscriptional Nuclear Regulation 

The formation of a pre-miRNA involves several factors, includ- 
ing phosphoprotein p53, estrogen receptor alpha (ERa), 
breast cancer 1 (BRCA1), protein p68, protein p72, splicing 
factor 2 (SF2/ASF), heterogeneous nuclear ribonucleoproteins, 
and KH-type splicing regulatory protein (Michlewski et al. 
2008; Wu et al. 2010; Suzuki and Miyazono 2011; Kawai 
and Amano 2012; Sundaram et al. 2013). These factors en- 
hance or inhibit miRNA maturation (fig. 4a). For instance, p53 
enhances Drosha activity to produce a miR-34 precursor 
(Tarasov et al. 2007), but it represses the maturation of miR- 
17-92 cluster (Yan et al. 2009). The binding of KH-type splic- 
ing regulatory protein to primary transcripts is essential for the 
processing of mir-198 (Sundaram et al. 2013). The docking 
of ERa to Drosha represses pri-miRNA processing of 
miR-1 6, miR-1 45, and miR-1 95 (Yamagata et al. 2009). 
Heterogeneous nuclear ribonucleoprotein A1 binds to the 
loop region of pri-miR-1 8a to facilitate slicing (Michlewski 
et al. 2008). Also, tissue-specific inhibition of pri-miR-7 pro- 
cessing by musashi RNA-binding protein 2 (MSI2) and human 
antigen R (HuR) proteins is reported in mammals (Choudhury 
et al. 2013). Some of these regulatory proteins act at multiple 
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Fig. 4. — as- and frans-acting regulatory elements during miRNA biogenesis (a) transcriptional regulation, where the genomic location of a miRNA 
determines its regulation. Green and brown bars indicate upstream promoter elements and repetitive sequences, respectively, (b) Posttranscriptional 
regulation: (i) Nuclear regulation, where the processing of pri-miRNA to pre-miRNA is determined by direct and indirect interactions with proteins, up- or 
downstream sequence elements and other factors (left), as well as the competition between microprocessors and spliceosome for a primary transcript that 
contains a segment of pre-miRNA at its exon-intron junction (right), (ii) Cytosolic regulation, where different factors affect pre-miRNA maturation (left) and 
degradation (right). See the text for the details. Hammer-headed red lines depict repression, arrow-headed green lines stand for promotion. Gray-dotted lines 
indicate interaction or alternative pathways. 



levels of miRNA biosynthesis pathway (interested readers are 
referred to Winter et al. 2009; Rinn and Huarte 201 1; Suzuki 
and Miyazono 2011; Finnegan and Pasquinelli 2013). The 
competition between miRNA processors and spliceosome 
has been shown in modulating the expression level of 
mature miRNAs (Ramalingam et al. 2014). Similarly, 52 
miRNAs, such as hsa-miR-202, hsa-miR-365, and hsa-miR- 
412 have active splice sites within a pri-miRNA, which can 
be a regulatory factor for their expression in tissue- and de- 
velopment-specific manner (Melamed et al. 201 3). Also, alter- 
native splicing may uncouple the expression pattern of 
clustered miRNAs from each other (Ramalingam et al. 2014). 

Autoregulation 

The structure of RNA itself has an important impact on a 
miRNA biogenesis pathway. miRNAs can regulate own bio- 
genesis through their secondary structure. For example, a 
miRNA cluster can promote or hinder miRNA processing 
through accessibility for miRNA processing machineries 
(Yang et al. 2009). Also, miRNAs are regulated by a feedback 
loop mechanism with their target (Yang et al. 2009), such as 
miR-57 and nob-1, or miR-7 and SF2/ASF (Wu et al. 2010; 
Zhao et al. 2010). In other example, let-7 binds at the 3'-end 
of its primary transcript and enhances its processing 



(Zisoulis et al. 2012). In addition, Drosha and Dgcr8 regulate 
each other (Han et al. 2009). 

Shuttle Regulation 

During the pre-miRNA transportation from the nucleus 
to cytoplasm, Exportin-5 protein protects pre-miRNA from 
degradation (Winter et al. 2009). The competition between 
pre-miRNAs and Dicer mRNA for Exportin-5 regulates the ma- 
ture miRNA homeostasis (Bennasser et al. 2011). Disturbed 
transport of a pre-miRNA to cytoplasm was observed in sev- 
eral cancer cell lines (Lee et al. 2008; Melo et al. 2010) and 
during viral infection (Bennasser et al. 201 1). 

Cytosolic Regulation 

miRNA maturation process can be further inhibited in cyto- 
plasm (fig. 4b); for example, Lin-28 binds to the terminal loop 
of let-7, and the 3'-end is polyuridylated by terminal uridyl 
transferases (TUT4/Zcchc1 1), thereby blocking Dicer process- 
ing (Heo et al. 2009). Another RNA-binding protein, MCPIP1, 
counteracts Dicer processing via cleavage of the terminal loop 
of a pre-miRNA (Suzuki et al. 201 1). Furthermore, phosphor- 
ylation of trans-activator RNA-binding protein, which is 
mediated by MAPK, enhances the stability of the miRNA- 
generating complex and results in an increase in miRNA 
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production and miRNA-mediated target silencing (Paroo 
et al. 2009). Loading miRNAs to Ago protein is a critical 
regulatory step (Suzuki and Miyazono 2011; Treiber et al. 
2012), which also determines the target specificity. Some 
studies showed developmental stage-specific abundance of 
one of the strands (either guide or passenger one) during 
the embryonic development of teleosts (Soares et al. 2009; 
Bizuayehu et al. 2012a). Recent studies show that miRISC 
can be bound by other transcripts such as circular RNAs 
(Hansen et al. 2013; Memczak et al. 2013), indicating 
additional level of regulation of mature miRNAs or different 
mechanisms of action. 

miRNA Decay 

Half-life of miRNA varies among tissues and miRNA types. It 
depends on target complementarity, 3'-end modification, cel- 
lular condition, and extracellular signaling (Katoh et al. 2009; 
Ruegger and GroBhans 2012). Analysis of miRNA turnover in 
mammalian embryonic fibroblasts showed the average half- 
life of miRNA was around 5 days in Dicerl -ablated cells 
(Gantier et al. 2011), which is by far greater than the half- 
life of mRNAs, 7.1 h on average (Sharova et al. 2009). 
However, other research showed faster decay of miRNAs; 
for example, miR-16 family stability is regulated in a cell 
cycle-dependent manner, in which the stability increases 
during the cell cycle exit and decreases during the re-entry 
(Rissland et al. 2011). 

miRNA Modifications and IsomiRs 

miRNA modifications are not stochastic. There are two types 
of miRNA modifications: Nucleotide modification and nucle- 
otide addition. The former one is an epigenetic mechanism; 
for example, adenosine (A) deamination results in conver- 
sion to inosine (I). A-to-l editing is the most common mod- 
ification; inosine has similar properties as guanosine (G) in 
base pairing, thus it can alter the pre-miRNA structure and 
mature sequences by affecting both miRNA processing 
and target recognition (Kawahara et al. 2007, 2008). For 
instance, A-to-l editing inside the seed sequence has been 
observed in some miRNAs, such as miR-151, miR-376a, miR- 
376b, and miR-368 (Kawahara et al. 2007; Garcia-Lopez 
et al. 2013). At least 6% of pri-miRNAs in mammals have 
A-to-l editing sites. The edited miR-376a was expressed 
in specific tissues, and it regulated target genes different 
than the unedited miR-376a (Kawahara et al. 2007). The 
edited precursors were removed at early postzygotic stages 
during mouse preimplantation development (Garcia-Lopez 
et al. 2013). 

Apart from A-to-l editing, mature miRNAs can be modified 
at the 3'-end through uridylation or adenylation (Katoh et al. 
2009; Chiang et al. 2010). This type of editing has been found 
in many miRNAs (Luciano et al. 2004; Kawahara et al. 2007, 
2008). It has important implications in miRNA biogenesis and 



target diversification, because it can affect the secondary 
structure of pri-miRNAs or miRNA:mRNA base pairing. 
Thus, it can create differential accumulation of mature 
miRNAs and target discrimination. The purpose of these 
modifications is context dependent; for example, adenyla- 
tion is required for selective stabilization of miR-122 in 
mouse liver (Katoh et al. 2009), whereas in THP-1 cell 
line, adenylation reduces effectiveness of miR-26a, miR- 
27a, and miR-122 (Burroughs et al. 2010). In vitro editing 
of pri-mir-142 in two positions remarkably reduced pre-mir- 
142 synthesis (Yang et al. 2006). This indicates that miRNA 
editing is one of the mechanisms that increase the reper- 
toire of miRNAs and their targets. 

Several mature miRNAs have size variants, termed as 
isomiRs. They are present in divergent species (Lee et al. 
2010; Li et al. 2011; Bizuayehu et al. 2012a; Humphreys 
et al. 2012; Wei et al. 2012; Yi et al. 2013). The origin of 
isomiRs is still not fully understood. Different mechanisms 
have been proposed, including degradation or imprecise 
cleavage of pre-miRNA during processing. However, the oc- 
currence of isomiRs is likely nonrandom. Although the random 
degradation of mature miRNAs by nucleases cannot be ex- 
cluded, the differential expression patterns and the observed 
target differences (Bizuayehu et al. 2012a; Humphreys et al. 
201 2; Wei et al. 201 2) suggest that biosynthesis of isomiRs is a 
regulated process. This hypothesis is supported by convergent 
results obtained in divergent species (Fernandez-Valverde 
et al. 2010; Cloonan et al. 2011; Wei et al. 2012). For exam- 
ple, nontemplate nucleotide additions at 3'-end, mostly A and 
U, but also C and G, have been reported in zebrafish (Wei 
et al. 2012), Atlantic halibut (Bizuayehu et al. 2012a), and 
blunt snout bream (Yi et al. 2013). A and U additions can 
stabilize or degrade mature miRNAs (Burroughs et al. 2010). 
Posttranscriptional miRNA modifications, resulting in a 
nontemplate nucleotide addition, involve a number of en- 
zymes including MTPAP, PAPD4, PAPD5, ZCCHC6, 
ZCCHC1 1, and TUT1 (Wyman et al. 201 1). 

Modifications occur not only at the 3'-end of a mature 
miRNA but also at its 5'-end. The majority of the 5' -end nu- 
cleotide alterations in Atlantic halibut miRNA were isomiRs 
(Bizuayehu et al. 2012a). It has been suggested that 5'-end 
size variations could result from the presence of multiple loci 
with different pre-miRNA structures (Starega-Roslan et al. 
2011). In the light of canonical miRNA:mRNA interaction 
and seed complementarity, the addition or truncation of nu- 
cleotides at the 5'-end can alter the target specificity of a 
miRNA. This has been shown in 5'-isomiRs of miR-101 
(Llorens et al. 2013) and miR-133a (Humphreys et al. 2012). 
However, Cloonan et al. (201 1) have shown that both canon- 
ical miRNAs and their isomiRs cooperate and have consider- 
able mRNA target overlaps. More studies are needed to 
uncover whether isomiRs are random degradation or rather 
regulated biosynthesis products, and whether the action of 
isomiRs is divergent or redundant. 



1916 Genome Biol. Evol. 6(8): 191 1-1937. doi:10.1093/gbe/evu151 Advance Access publication July 22, 2014 



miRNA in Teleost Fish 



GBE 



miRNA Target Site Polymorphism and Posttranscriptional 
Modifications 

Polymorphism in DNA sequence of the target site influences 
miRNA-target interaction through stabilizing or destabilizing 
the existing miRNA target sites, or creating new target sites 
(Georges et al. 2007). miRNA target site polymorphism can 
result in phenotypic variation and disease conditions (Saunders 
et al. 2007; Ziebarth et al. 2012). For instance, a single-nucle- 
otide polymorphism (SNP) at the 3'-UTR of mystotin gene in 
Texel sheep allows binding by miR-1 and miR-206, which in 
effect creates muscular hypertrophy (Clop et al. 2006). 

In teleosts, target site polymorphism is found in a number 
of processes, such as muscle development and regeneration, 
photoreceptor morphogenesis, immune response, or cranio- 
facial development, and it can have functional effects (Loh 
et al. 201 1 ). In the three-spined stickleback (Gasterosteus acu- 
leatus), polymorphic target sites were found at the 3'-UTR of 
Glyceraldehyde-3-phosphate dehydrogenase, which has two 
alleles. One allele was predicted as a target of miR-2888, miR- 
705, and miR-2305, whereas the other allele was targeted by 
miR-1 777b (Chaturvedi et al. 2014). In Lake Malawi cichlids, 
SNP density in the predicted miRNA target sites was higher 
than in the flanking regions. Allele frequency analysis and lin- 
eage specificity of these sites suggested contribution of 
miRNA target site polymorphism to species diversification 
(Loh et al. 201 1). Further identification of miRNA target site 
polymorphism in other teleosts will help to understand the 
role of miRNA in teleost evolution. 

RNA editing has an effect not only on miRNA processing 
and target diversification; A-to-l and C-to-U editing can also 
disrupt legitimate miRNA target sites. It has been found that 
RNA editing sites are highly enriched at "seeds" of miRNA 
target sites; and this can create new miRNA target positions or 
disrupt the existing ones (Gu et al. 2012). Editing can alter 
RNA secondary structure, which can affect the accessibility of 
miRNA-binding sites (Brodersen and Voinnet 2009). 

Evolutionary Constraints and Teleost Specificity of 
miRNAs 

Selection pressure on mature miRNA is immense. How the 
potency of miRNA fine tuning of diverse biological pathways 
is rendered in the course of evolution? This capacity is shaped 
by various selective forces (internal and external), such as pre- 
cursor structure, RNA-binding proteins, target-binding sites, 
target selection, and decoys. The "minimal sufficient" struc- 
tural requirement of miRNA are as follows: 1) A pre-miRNA 
must form a stem structure homoduplex, meaning that two 
segments have to show complementarity with low free 
energy; 2) an miRNA must fulfill the minimum requirement 
to be loaded to Argonaute protein, and in the case of animals, 
no perfect complementarity to Ago catalytic center is re- 
quired; and 3) miRNA's seed sequence must have binding 
site(s) on its target(s). Moreover, miRNA evolution is affected 



by decoys, such as circular RNAs, long noncoding RNAs, and 
other small RNAs (Hansen et al. 2013; Memczak et al. 2013). 
In addition, miRNA can be bound by competing endogenous 
RNAs (ceRNAs). Therefore, miRNA loci are under both positive 
and negative selection pressure (Kosik 2013). 

Numerous miRNAs are common among divergent animal 
species (Niwa and Slack 2007). Several lines of evidence indi- 
cate that the speciation of metazoans has been accompanied 
by emerging novel miRNAs (Hertel et al. 2006; Tarver et al. 
201 3) and that the majority of the inventions have occurred in 
vertebrates (Heimberg et al. 2008). Teleosts have additional 
copies of miRNAs as a result of duplicate retention following 
the teleost-specific whole-genome duplication and gene du- 
plication events. A study on localization of selected pri-miRNA 
duplicates in zebrafish has been performed (He et al. 201 1), 
but no information is available on the regulation and func- 
tional characterization of miRNA duplicates, that is mature 
miRNA sequences originating from different loci. However, 
it is possible that cell- or tissue-specific transcriptional factors 
may induce differential expression of miRNA duplicates. For 
instance, human miR-365 has two copies, one with an active 
splice site, whereas the other one is intronic; this suggests 
regulator modulation of the expression of paralogous 
miRNAs (Melamed et al. 2013). Previous studies on protein- 
coding genes demonstrated that several paralogs, which 
evolved as a result of teleost-specific whole-genome duplica- 
tion, had distinct pattern of expression and attained subfunc- 
tionalization, neofunctionalization, or loss of their function in 
the course of evolution (Brunet et al. 2006; Crow et al. 2006). 
Furthermore, some groups of teleosts underwent genome 
reduction, such as members of orders Tetraodontiformes 
and Pleuronectiformes (Brainerd et al. 2001; Venkatesh 
2003); thus, studying miRNAs in these species in comparison 
to species that retain duplicates can expand our understand- 
ing of the genome duplication event. Functional characteriza- 
tion of miRNA duplicates is important in this context. 

miRNA in Teleost Development 

Early Development 

miRNAs are involved in regulation of early developmental tran- 
sitions. A number of miRNAs has temporally defined expres- 
sion patterns, such as those expressed during maternal- 
to-zygotic transition (MZT) and metamorphosis (table 2). 
MZT is a two-step process of removal of a subset of maternal 
mRNA and proteins followed by the initiation of zygotic mRNA 
transcription (Tadros and Lipshitz 2009). It is regulated by the 
two networks, maternal and zygotic. Maternal mRNA stability, 
translation, and localization are the three features that deter- 
mine control of early embryogenesis. miRNA expression pro- 
filing experiments have shown the probable role of miRNA in 
destabilization of maternal transcripts (Wienholds et al. 2005; 
Tani et al. 2010; Bizuayehu et al. 2012a). mRNA degradation 
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pathway through the miRNA mediation occurs after the MZT 
(Barckmann and Simonelig 2013). miR-430, which is highly 
expressed during the blastula stage, is involved in maternal 
transcript clearance (Giraldez et al. 2006). However, other 
miRNAs such as miR-34, miR-200a, miR-200b, and miR-206 
are also found abundantly during the embryogenesis. For ex- 
ample, miR-206 is expressed both maternally and zygotically 
in zebrafish, and it is essential in controlling cell movements 
during the gastrulation (Liu et al. 2012). 

Several studies demonstrated the role of miRNAs in the 
progression of teleost embryonic development. Zebrafish em- 
bryos and larvae lacking zygotic Dicerl had slow growth rate 
and survived only for 2 weeks (Wienholds et al. 2003). 
Similarly, maternal and zygotic Dicer mutant zebrafish em- 
bryos had morphogenetic defects during the gastrulation, 
brain formation, somitogenesis, and heart development 
(Giraldez et al. 2005). In the latter study, injections of 
mature miR-430 into the Dicer-deficient embryos partially res- 
cued the gastrulation and reduced brain ventricle morphogen- 
esis defects, indicating its role in the processes. In addition, 
miR-430 is transcribed after the zygotic genome activation by 
maternally stocked transcriptional factors, such as Nanog, 
Pou5f1, and SoxB1 (Lee et al. 2013), indicating zygotic 
origin of this miRNA. However, miR-34 is a maternal miRNA 
involved in early neural system development (Soni et al. 2013). 
The roles of other maternally stocked miRNAs, such as miR-24, 
miR-30, miR-126, miR-146, and miR-221 (Ma et al. 2012; 
Juanchich et al. 2013) remain to be uncovered. 

Organogenesis 

Several species of miRNAs have been characterized during 
teleost organogenesis (table 2). 

Asa part of the transcript pool, miRNAs create a context for 
the organ to be formed. Rudiment formation of any organ 
needs organized causality, which starts with signaling fol- 
lowed by consequential changes in transcripts pool manage- 
ment. These sequential actions require superseding the 
transcripts, buffering the noise from unintended transcripts, 
and shaping the transcriptional output to fit the context or 
keep homeostasis. miRNAs help to establish discrete domains 
of gene expression during organogenesis. Below we discuss 
known or anticipated roles of some miRNAs in the formation 
of different teleost tissues. 

Brain Formation 

Diverse types of miRNAs are present in distinct regions of brain 
(Kapsimali et al. 2007), implying constricted function in a 
given region. In teleosts, conserved brain-specific miRNAs 
are found in divergent species (Soares et al. 2009; Xia et al. 
2011; Bizuayehu et al. 2012b; Zhu et al. 201 2; Xu et al. 201 3). 
Zebrafish maternal and zygotic Dicer mutant embryos showed 
proper developmental progression of neural plate to neural 
rod; however, a considerable impairment in neural 



development was observed in the formation of the neurocoel 
and neural tube, as well as reduction of the brain ventricles 
and lack of distinct brain regions were found. These defects 
were partially rescued by injection of a preprocessed miR-430 
family to the mutant (Giraldez et al. 2005). 

The expression of brain miRNAs depends on the cell status; 
for example, miR-92b is widely expressed in proliferative 
neural cells regardless of the fate of these cells, whereas 
miR-124 is expressed in differentiated neurons only. In con- 
trast, miR-9 and miR-135c are expressed in both cell types 
(Kapsimali et al. 2007). In that study, the authors also 
showed that miR-181a and b were expressed specifically in 
retina cells. 

miRNAs have brain-organizing activity; for instance, miR-9 
is expressed selectively in late embryonic neural tube by spar- 
ing the midhind brain to define the boundary (Leucht et al. 
2008). Other study showed that loss of miR-7 could result in 
specific reduction of midbrain size without affecting the tel- 
encephalon at the anterior tip of the brain (Memczak et al. 
2013). In summary, localized, transient, and constitutive ex- 
pression of miRNAs in teleost brain indicates their function in 
brain morphogenesis and maintenance of distinct subregions 
and cell types. 

Eye Formation 

The embryonic origin of teleost eye is similar to other verte- 
brates; however, vision in teleost depends on ecological niche 
and behavior of a species. The vision procures adaptation to 
the environment, thus eye tissue-specific gene expression 
guides this adaptation. Eye development has been well char- 
acterized in zebrafish using morphology, gene expression, and 
in situ labeling (Fadool and Dowling 2008; Gestri et al. 2012). 
During gastrulation, the middle part of anterior neural plate is 
destined to be an eye field. This field is under the control of 
different signaling pathways, which influence the develop- 
ment of the eye. Among them, Wnt signaling pathway de- 
fines regions of anterior neural plate including eye field and 
migration of eye field cells, and promotes eye formation 
(Cavodeassi et al. 2005). The induction of the eye field is fol- 
lowed subsequentially by a formation of optic cup through 
invagination of optical vesicles, change of the optic stalk to 
optic nerve and retina, and by the closure of choroid fissure 
(Gestri et al. 2012). This remodeling is regulated by many 
transcriptional factors, such as Meis2, Mitf, Pax2, Pax6, 
Six3a, Vax1, and Vax2, and signaling pathways, such as Fgf, 
Hh, Shh, and Wnt (Macdonald et al. 1997; Chow and Lang 
2001; Cavodeassi et al. 2005; Conte et al. 2010). 

Several miRNAs, such as miR-96, miR-124a, miR-181a, 
miR-181b, miR-182, miR-183, miR-184, and miR-2 04 are ex- 
pressed in eye of zebrafish embryo (Cavodeassi et al. 2005). 
Similarly, a number of miRNA species was identified in Asian 
seabass (Lates calcarifer) eye (Xia et al. 201 1). Spatial localiza- 
tion of miRNAs revealed cell type- and developmental stage- 
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specific expression patterns (Kapsimali et al. 2007). For in- 
stance, miR-181a and b were expressed specifically in retina 
cells (Kapsimali et al. 2007). The authors also showed re- 
stricted expression of miR-92b and let-7b in the ciliary mar- 
ginal zone of the retina and complete absence of these 
miRNAs in mature retinal neurons. miR-9 was expressed in 
mature amacrine cells of the inner nuclear layer and in matur- 
ing cells of ciliary marginal zone of the retina. miR-30a, miR- 
184, and mir-204 were localized in lens. Conte et al. (2010) 
have showed that miR-204 targets meis2 and modulates Pax6 
transcriptional pathway in medaka Oryzias latipes. Using mor- 
pholino-based knockdown approach, they demonstrated that 
the depletion of miR-204 resulted in a number of eye devel- 
opment malformations including eye cup impairment, small 
eyed embryos, impaired lens development, defect in lens ep- 
ithelial cells patterning, misplacement and disorganization of 
primary fiber cells, lens herniation, and failure of optic fissure 
closure. The evolutionary conservation of Pax6-miR-204 path- 
way is demonstrated in mouse ocular tissues (Shaham et al. 
2013). Other miRNAs are found in keeping homeostasis of 
cells of the eye; for example, let-7 maintains Muller glia cells 
in a differentiated state (Ramachandran et al. 2010). 

Given the need for the precise and intricate regulation of 
gene expression during eye development and the adaptive 
significance of eye, further work on developmental and phys- 
iological functions of miRNAs in the eye of various teleosts 
from different ecological niches would disclose the adaptive 
role of miRNAs in teleosts by using visual system as a model. 

Muscle Formation 

Muscle formation in teleosts begins during the embryonic 
development with the formation of precursor myogenic 
cells. These myogenic cells differentiate into a myotome, 
which has four cell lineages: Muscle pioneers, slow muscle, 
fast muscle, and medial fast muscle (Johnston 2006). 
Specification of these lineages is regulated by myogenic reg- 
ulatory factors (MRFs). MyoD, a member of MRF, controls 
several downstream genes involved in myogenesis. MyoD is 
directly regulated by miR-203b in Nile tilapia (Oreochromis 
niloticus) (Yan, Guo, et al. 2013). In zebrafish, myotube pro- 
duction ceases at 40% of the total body length and the tran- 
sition from hyperplasia to hypotrophy is facilitated by miRNAs, 
including let-7, miR-19, and miR-130 families (Johnston et al. 
2009). 

During skeletal myogenesis, myoblasts differentiate into 
slow-twitch or fast-twitch muscle fibers. These lineage-specific 
pathways are established by the activity of either slow-twitch 
restricted or fast-twitch restricted genes, as well as transcrip- 
tion factors (Chauvigne et al. 2005; Elworthy et al. 2008). 
Lineage-restricted expression of miR-499 leads to the estab- 
lishment and maintenance of slow-twitch muscle fibers 
through repression of Sox6, which promotes fast-twitch 
muscle differentiation; this mechanism is conserved among 
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vertebrates (Wang et al. 2011). Several miRNAs have been 
identified in skeletal muscle of common carp, and some of 
them are unique for teleosts (Yan et al. 2012). In zebrafish 
embryos, miR-1 and miR-133 were implicated in shaping sar- 
comere actin organization (Mishima et al. 2009). Further ex- 
ploration of the function of these unique miRNAs is important 
to understand muscle formation in teleosts. The role of miRNA 
in muscle growth is addressed in another section. 

Cardiovascular Formation 

At the early stage of teleost organogenesis, beating heart is a 
recognized developmental stage because of its visibility ahead 
of other discernible organs. The heart differentiates from two 
distinct cardiac progenitor cells termed the first heart field and 
the second heart field (De Pater et al. 2009; Grimes et al. 
2010). miRNAs participate in the regulation of migration of 
the heart fields to the midline; for example, miR-218a-1/2 
titrates roundabout homolog 1 (robol) to regulate endocar- 
dial migration through vascular endothelial growth factor 
(Vegf) signaling (Fish et al. 201 1). In zebrafish, the first heart 
field gives rise to heart tube and progenitors of the second 
heart field differentiate to form smooth muscle and myocar- 
dium (Hami et al. 201 1). The heart tube loops, tightens, and 
forms atrium and ventricle, two chambers delineated morpho- 
logically, molecularly, and functionally. miR-1 38 knockdown 
during zebrafish cardiac development resulted in defects in 
elongation of ventricular cardiomyocytes and in early cardiac 
looping (Morton et al. 2008). In the further analysis, the 
authors showed that miR-1 38 had a restricted expression in 
atrioventricular canal, and it targeted atrioventricular canal 
domain genes by regulating retinoic acid synthesis and 
direct repression of chondroitin sulfate proteoglycan 2 
(cspg2). In zebrafish during heart development, miR-218a in- 
teracts with transcriptional factor Tbx5 (Chiavacci et al. 2012), 
which is necessary for endocardial cell differentiation and 
valve tissue formation (Camarata et al. 2010). Also, miR-21 
is crucial in regulation of heart valve formation by modulating 
the expression of sprout homolog (sprout), programmed cell 
death 4 (pdcd4K), and phosphatase and tensin homolog B 
(ptenb) (Banjo et al. 2013). miR-1 26 knockdown resulted in 
collapsing lumens and compromised endothelial tube organi- 
zation by repressing sprouty-related EVH1 domain-containing 
protein 1 (spredl) and phosphoinositide-3-kinase, regulatory 
subunit 2 (pik3R2) to promote Vegf signaling during zebrafish 
vascular development (Fish et al. 2008). Similarly, Lalwani et al. 
(2012) using knockdown and overexpression assays showed 
that miR-1 42 -3p targeted cdh5 and influenced vascular integ- 
rity, remodeling, and angiogenesis. Also, let-7 family, miR- 
20b, miR-31, miR-221, and miR-1 81 a promote angiogenesis 
and lymphangiogenesis in zebrafish (Biyashev et al. 2012; 
Nicoli et al. 2012; Dunworth et al. 2013). Studies in mammals 
documented several interesting functions of miRNAs during 
cardiac morphogenesis; similar studies in teleosts would be 



valuable to determine miRNA functional conservation in 
heart development. 

Gametogenesis 

Primordial germ cells (PGCs) are the carrier of the genetic 
information from one generation to another; therefore, mo- 
lecular events during their developmental progression must be 
strictly regulated to ensure a stable transmission of genetic 
information to future generations. In model fishes, gameto- 
genesis starts from asymmetric mitotic divisions of PGCs, 
which are specified very early during embryogenesis. PGCs 
migrate to the future genital ridges and become gonocytes, 
then during sex differentiation, they transform to spermato- 
gonia or oogonia (Lubzens et al. 2010; Schulz et al. 2010). 

PGC specification necessitates the suppression of somatic 
lineage programs. Depending on the species, PGCs are tran- 
scriptionally inert at the beginning (Nakamura and Seydoux 
2008; Venkatarama etal. 2010). PGCs have differential stock- 
pile of proteins and transcripts compared with that of somatic 
cells. One of the mechanisms securing selective mRNA tran- 
script profile in PGCs has been discovered in zebrafish 
(Mishima et al. 2006). miR-430, a major "clearance" miRNA 
during the early embryonic development, suppresses some 
transcripts, such as nanos, tdrd7, and hub in somatic cells 
but not in PGCs (Mishima et al. 2006; Kedde et al. 2007; 
Mickoleit et al. 201 1). Nanos, tdrd7, and hub are necessary 
for proper migration, maintenance, and survival of PGCs 
(Koprunner et al. 2001; Mickoleit et al. 2011), and Dead 
end (Dnd), an RNA-binding protein, protects 3'-UTR-binding 
sites of these transcripts from miR-430-mediated repression in 
PGCs (Kedde et al. 2007). Morpholino-mediated knockdown 
of dnd in several teleost species leads to removal of protection 
of key PGC-specific transcripts from miR-430-guided suppres- 
sion, and consequently PGC development is arrested, germ- 
line lineage is lost, and the developing individuals are 
irreversibly sterile (Weidinger etal. 2003; Fujimotoetal. 2010). 

After specification in an early embryo, PGCs migrate to the 
future gonadal ridge. This migration is guided by a chemokine 
Sdf-1 signaling from the neighboring somatic cells recognized 
by CXCR4 receptor (Doitsidou et al. 2002; Knaut et al. 2003). 
This mechanism, found for the first time in zebrafish 
(Doitsidou et al. 2002), is well conserved across the investi- 
gated vertebrates (Stebler et al. 2004). In this migratory route, 
the role of miR-430 has been demonstrated in the clearance 
of sdf-la mRNA from previous expressing domains to ensure 
correct migration of PGCs (Staton et al. 201 1). 

miRNAs are essential for proliferation and maintenance of 
germ cell-supporting somatic cells, such as Sertoli and Leydig 
cells in the testis, and follicle cells in the ovary. Sertoli cell 
number is a limiting factor in sperm production in a sexually 
mature fish (Schulz et al. 2010). In Sertoli cell-specific Dicer 
conditional knockout mouse model, miR-1 25a-3p, miR-872, 
and miR-24 have role in translational control during 
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spermatogenesis (Papaioannou et al. 2011). Rakoczy et al. 
(2013) reported up to 2.6-fold increase in the number of 
Leydig cells in miR-140-3p ablated mice. Also, miR-202-5p/ 
3p transcripts were identified as potential regulators of 
mouse embryonic gonad differentiation with strong expres- 
sion in Sertoli cells (Wainwright et al. 2013). Multiple miRNAs 
are implicated in granulosa cells apoptosis, steroidogenesis, 
and cell proliferation (Donadeu et al. 2012; Hu et al. 2012; 
Yang et al. 2012; Yin et al. 2012; Zhang et al. 2013b). 
Functional involvement of miR-224, miR-378, and miR-382 
in regulation of aromatase expression during follicle develop- 
ment has been demonstrated, and miR-21 promoted follicular 
cell survival during the ovulation (Donadeu et al. 2012). 
No such data exist in fish despite the fact that in seasonally 
reproducing teleost species, proliferation and apoptotic pro- 
cesses in gonads are orchestrated by a complex regulatory 
network (Nagahama 1994; Chaves-Pozo et al. 2005; 
Almeida et al. 2008). 

Skeletogenesis 

Skeletal formation includes cartilaginous state (chondrifica- 
tion) and ossified state (ossification) and is regulated by several 
ubiquitous and specific genes. The major part of skeleton orig- 
inates from the neural crest, lateral plate mesoderm, paraxial 
mesoderm, and notochord. Mesenchymal cells become chon- 
drocytes, osteoblasts, or other skeletal cells depending on 
transcriptional factors and signaling pathways (Karsenty and 
Wagner 2002). Transcriptional factor sox9, for example, 
is essential for morphogenesis of condensation and cartilage 
differentiation in zebrafish (Yan et al. 2002) and regulates the 
expression of miR-140 (Nakamura et al. 2012). However, miR- 
140 acts independently of sox9 in the regulation of palatal 
skeleton development by modulating pdgf-receptor alpha, 
which is required for migration of palatal precursor and 
neural crest cells (Eberhart et al. 2008). Similarly, Runx2 is a 
transcriptional factor that regulates osteoblast differentiation 
(Flores et al. 2006). Huang et al. (2010) showed direct nega- 
tive regulation of Runx2 by miR-204/miR-21 1 in stroma and 
myoblast cell lines. These results indicate the importance of 
miRNAs during bone and cartilage formation but still those 
few miRNAs that are profiled as skeleton specific are not char- 
acterized functionally, particularly in the context of the pre- 
vailing problem of skeletal deformities experienced in fish 
production (Silverstone and Hammell 2002; Bardon et al. 
2009). 

miRNA in Growth and Regeneration 
of Teleosts 

Growth 

Muscle cell-type specification requires Hedgehog, fibroblast 
growth factor, and retinoic acid signaling pathways, as well 
as T-box genes (Lewis etal. 1999; Ochi and Westerfield 2007). 



Both coordinated spatiotemporal action of regulatory factors 
and the removal of certain domains from the previously 
expressing cells are essential in the cell specification. In zebra- 
fish, miR-21 4 enhances the cellular response to Hedgehog 
signaling and facilitates strict specification of muscle cell 
types through the repression of su(fu), which is a negative 
regulator of Hedgehog signaling (Flynt et al. 2007). Bone mor- 
phogenetic protein and transforming growth factor beta 
(TGFp) signaling pathways modulate the activity of some 
proteins, such as Smadl, p68, Drosha, or Dgcr8, which are 
involved in miRNA biogenesis (Suzuki and Miyazono 2011). 
Once the cell commitment occurs, differential regulation of 
various factors promotes the differentiation and proliferation, 
which allow subsequent tissue development. 

Fish growth is modulated by spatiotemporal expression of 
various genes. In rainbow trout (Oncorhynchus mykiss) alevins, 
the transition from endogenous (yolk sac stage) to exogenous 
feeding changes the expression of metabolic genes and 
miRNAs (Mennigen et al. 2013), indicating the participation 
of miRNAs in metabolic pathways of teleost fish during their 
early growth. Huang et al. (2011) have reported a negative 
feedback circuit in which insulin-like growth factor 1 (IGF-1) 
promotes miR-133 expression, which, in turn, represses IGF-1 
receptor (IGF-1 R) affecting skeletal myogenesis. In Nile tilapia 
(O. niloticus), miR-206 targets IGF-1 and inhibits its action 
(Yan, Zhu, et al. 2013), indicating the importance of 
miRNAs in hypothalamic-pituitary pathway. 

Growth of skeletal muscle mass in fish occurs through hy- 
pertrophy and hyperplasia of muscle fibers, which absorb 
myoblasts differentiated from myogenic precursor cells 
(Johnston 2006). Transcriptional factors, signaling proteins, 
and ncRNAs are the determinants of muscle mass formation 
from the first cell commitment to the last fusion stages. These 
steps include myogenic progenitor cells specification, activa- 
tion, proliferation, cell cycle exit, differentiation, migration, 
and fusion. In early embryogenesis of teleosts, specifically 
during gastrulation, the commitment of somatic cells to be 
myogenic cells is ignited by MRFs, which also stimulate 
muscle-specific miRNA biogenesis (Sweetman et al. 2008). 
On the other hand, miRNAs regulate the level of MRFs. 
In zebrafish, Goljanek-Whysall et al. (2011) demonstrated 
elevated expression of pax3 in dermomyotomal progenitors, 
but its downregulation in the MRF-expressing myotome. In 
their model, MRFs activate miR-1/miR-206 in the committed 
myoblasts, and these miRNAs target residual pax3 during the 
progenitor-to-myoblast transition, and they control transi- 
tional timing by repressing pax3. 

Several miRNAs regulate teleost skeletal muscle growth. 
For example, Huang et al. (2012) reported significant differ- 
ential expression of skeletal muscle miRNAs between fast- 
growing and slow-growing strains of Nile tilapia, indicating 
possible application of miRNAs as selection markers for aqua- 
culture industry. In another study, four miRNAs (miR-1, miR- 
27a, miR-1 33a, and miR-206) were differentially expressed 
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during skeletal muscle development of Nile tilapia (Yan et al. 
2012a). Similarly, by comparing skeletal muscle of different 
stages (larvae, 1-, and 2-year old) of common carp (Cyprinus 
carpio), Yan et al. (2012) reported an increase in miR-1, miR- 
21, miR-1 33a-3p, and miR-206 expression with age. Earlier in 
situ hybridization study in zebrafish showed localization of 
expression of these miRNAs in skeletal muscle (Wienholds 
et al. 2005). In adult zebrafish, Johnston et al. (2009) demon- 
strated significant differences in miRNA expression between 
two fast muscle phenotypes, myotube recruiting and ceased. 
Further investigation of the relation between miRNAs and 
MRFs is essential in understanding skeletal muscle growth. 

Growth depends on feed intake and is influenced by a 
number of intrinsic and extrinsic factors (Hoskins and 
Volkoff 2012), which modulate the expression of miRNAs 
and their targets. For example, in rainbow trout alevins, a 
change from endogenous to exogenous feed is associated 
with the high expression of miR-1 43 and inverse expression 
of its target abhd5 (Mennigen et al. 2013), a gene that acti- 
vates adipose triglyceride lipase. Feed intake depends on ap- 
petite, which is regulated by hypothalamus. Several appetite 
stimulators/orexigenes (e.g., orexins, neuropeptide Y, and 
ghrelin) and inhibitors/anorexigenes (e.g., cholecystokinin, 
leptin, and amylin) have been isolated in fish (Volkoff et al. 
2005). The roles of miRNAs in orexigenic and anorexigenic 
pathways would be the important area of investigation to 
utilize miRNAs as markers for selective breeding programs. 

Sexual dimorphism in growth is observed in many teleost. 
This difference is largely associated with reproduction-related 
features, such as timing of maturation, territorialism, orna- 
mentalism, courtship, nest guarding, and nursing, in which 
the amount of spent energy affects the growth of one of 
the sexes particularly (Hendry and Berg 1999). Various studies 
have shown the difference in expression of genes implicated 
in sexual growth dimorphism; however, few studies have 
been performed on sexually dimorphic expression of 
miRNAs in fish tissues (Bizuayehu et al. 2012b) not directly 
related to the growth of fish. Scientific knowledge of the ge- 
netic basis of growth traits provides an avenue for improve- 
ment in aquaculture, and miRNAs can be assistive in such the 
endeavor. 

miRNAs are implicated in nutrient metabolism. Fasting and 
re-feeding experiment in rainbow trout showed significant 
upregulation of miR-1 22 and miR-33 together with cpf7a 
and cptlb in liver, suggesting lipogenic role of miRNAs at 
multiple levels of the hepatic intermediary metabolism 
(Mennigen et al. 2012). Twenty-seven growth-related 
miRNAs were identified in blunt snout bream Megalobrama 
amblycephala, including miR-23b, miR-92, and miR-462, 
which were expressed abundantly in slow-growing groups 
compared with fast-growing groups (Yi et al. 2013). The pre- 
dicted targets of these miRNAs were involved in metabolic 
pathways. 



Regeneration 

Teleosts grow continuously throughout their life, and they are 
able to regenerate organs and appendages, including spinal 
cord, heart, retina, scales, and fins. Although similar molecular 
and cellular processes exist in both mammals and fish, regen- 
erative capacity of adult mammals is limited when compared 
with fish. Concentration or spatiotemporal availability of de- 
velopmental regulators might constitute a difference in regen- 
eration capability between mammals and teleosts (Kawakami 
et al. 2006). A number of studies have been conducted to 
elucidate biological mechanisms that govern regeneration, 
but still there are unidentified signals that initiate and regulate 
certain regenerative steps. Epimorphic regeneration requires 
signaling pathways, such as Wnt, fibroblast growth factor, 
retinoic acid and Hedgehog, or notch (Poss et al. 2003; 
Kawakami et al. 2006; Blum and Begemann 2012; 
Weidinger et al. 2013). The requirements for a signal can be 
similar or variable among the regenerative tissues (Duszynski 
et al. 2013). As modulators of transitional events, miRNAs 
participate in the regulation of regeneration process through 
temporal clearance of unneeded transcripts. 

Well-investigated regeneration of caudal fin in teleosts 
has three main discernible stages: 1) wound healing, closing 
of wound by the migration of epithelial cells; 2) blastema 
formation, disorganization of mesenchymal tissue followed 
by formation of a mass of undifferentiated, proliferating mes- 
enchymal progenitor-like cells; and 3) regenerative out- 
growth, proliferation, and differentiation of blastemal cells 
(Poss et al. 2003). Numerous miRNAs are down- or upregu- 
lated during the fin regeneration. For instance, miRNA micro- 
array experiment showed the upregulation or maintenance of 
6 miRNAs and downregulation of 1 6 miRNAs by Fgf signaling 
(Yin et al. 2008). Massive differences in miRNA expression are 
reported between intact, amputated, and regenerating fins of 
zebrafish (Thatcher et al. 2008). This report indicated that 
bmp3, hsp60, and msxb genes had role in regeneration and 
were predicted targets of up- or downregulated miR-200b, 
miR-2/miR-338, and miR-301, respectively. In gain- and loss- 
of-function experiments, Yin et al. (2008) showed that the 
regulated depletion of miR-1 33 resulted in effective fin regen- 
eration. Similarly, miR-203, which has binding sites on the 
3'-UTR of transcriptional factor lefl, has been significantly 
downregulated during fin regeneration (Thatcher et al. 
2008). Lefl is involved in blastema formation and marks the 
basal epidermal layer and distal blastema (Poss et al. 2000). 
miR-203 represses left! and in consequence blocks the fin 
regeneration (Thatcher et al. 2008). In contrast, loss of miR- 
203 results in abundance of lefl and fin overgrowth, indicat- 
ing the importance of miR-203 not only in blastema formation 
but also for proper termination of the regeneration process 
(Thatcher et al. 2008). These findings indicate the reprogram- 
ming of miRNA expression in a tissue in response to the 
regeneration program. 
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Cardiac tissue regeneration in zebrafish occurs through 
cardiomyocyte dedifferentiation, transdifferentiation, and pro- 
liferation (Jopling et al. 2010; Zhang et al. 2013c). Alteration 
of expression of several miRNAs during cardiac regeneration 
(upregulation of 10 and downregulation of 8 miRNAs) was 
observed at 7 days postamputation (Yin, Lepilina, etal. 2012). 
Further evaluation of miR-133 in this study indicated that miR- 
133 had several targets, among them mpsl and cx43, which 
are essential for the regeneration process. Given that miR-1 33 
has regulatory role in skeletal muscle proliferation (Chen et al. 
2006), the repression of miR-133 during heart regeneration 
may indicate reprogramming. 

miRNAs are implicated in neuronal regeneration, including 
central nervous system and retina. Repressive action of miR- 
133b toward rhoA mRNA, which inhibits axon regrowth, has 
been demonstrated during spinal cord regeneration (Yu et al. 
2011). In addition, the authors have demonstrated that 
miR-1 33b expression in medial longitudinal fascicle, superior 
reticular formation, and intermediate reticular formation neu- 
rons is essential for full locomotor recovery after a spinal cord 
injury. During retinal regeneration, dedifferentiation of Muller 
glia into a cycling population of progenitor cells enables the 
injured retina to restore its function (Ramachandran et al. 
201 0). This process is partially stimulated by asclla regulation 
of Lin-28 protein, which decreases let-7 level. Lin-28 and let-7 
have inverse regulatory link (Rybak et al. 2008). Thus, inhibi- 
tion of let-7 promotes the expression of genes that are nec- 
essary for retina regeneration in zebrafish (Ramachandran 
etal. 2010). 

The mechanisms of regeneration are complex and involve 
multiple signaling molecules, transcriptional factors, and 
genes. So far, few miRNA targets have been validated. 
Although many pathways have shown overarching conserva- 
tion among the regenerative tissues, further exploration 
of tissue specificity and conservation of miRNAs among the 
regenerative organs may provide better knowledge of regen- 
erative mechanisms. 

miRNA in Teleost Reproduction 

Oogenesis and Spermatogenesis 

Studies in mammals indicate that miRNAs are crucial in oo- 
genesis and spermatogenesis (Takada et al. 2009; Yadav and 
Kotaja 2014). The balance between self-renewal and differ- 
entiation of spermatogonia! cells is critical in seasonally repro- 
ducing fishes to define the start and the end of reproductive 
activity; thus, the maintenance of this balance requires post- 
transcriptional regulation of a number of genes. However, no 
information is available on miRNA functions in the spermato- 
gonia! phase of teleost spermatogenesis. In miRNA microarray 
experiment, differential expression of 1 3 miRNAs was found 
at previtellogenesis, vitellogenesis, late vitellogenesis, and 
maturation stages during oogenesis in rainbow trout. 



miRNA targets, important in oocyte maturation, growth, 
development, and maturational competence, that is the ability 
of oocyte to resume meiosis, were also predicted (Juanchich 
et al. 2013). However, this and other few studies on miRNA 
in teleost gonadal development (Bizuayehu et al. 2012b; 
Abramov et al. 201 3) give an insight into miRNA developmen- 
tal profile rather than decipher cell type-specific functions. 

Many circulating endocrine and locally acting paracrine and 
autocrine factors regulate oogenesis and spermatogenesis 
(Lubzens et al. 201 0; Schulz et al. 201 0). Estrogens are impor- 
tant in teleost spermatogonia! self-renewal and in the transi- 
tion of type A to type B spermatogonia (Schulz et al. 2010). 
There is a complex interplay between estrogens and miRNAs. 
Estrogens regulate transcription of some miRNAs, such as 
miR-21 and miR-221; in contrast, other miRNAs such as let- 
7, miR-22, miR-1 96b, or miR-206 target estrogen receptor 
alpha transcript (Cochrane et al. 2011). In mammals, miR- 
383 regulates 1 7p-estradiol release from granulosa cells (Yin 
et al. 2012). Also androgens regulate the expression of some 
miRNAs, such as miR-22, miR-1 22, and mir-125b (Cochrane 
et al. 2011); however, direct posttranscriptional regulation 
of androgens by miRNAs is not elucidated yet. The effect of 
masculinization treatment with either a synthetic androgen 
(17-a-methyl testosterone) or an inhibitor of cytochrome 
P450 aromatase (Fadrozole) on miRNA expression has been 
studied in Atlantic halibut; masculinization treatment resulted 
in differential expression of let-7a, miR-1 9b, miR-24, and 
miR-202-3p in gonads (Bizuayehu et al. 2012b). Future 
work is needed to investigate the role of miRNAs in hor- 
mone-secreting gonadal cells. 

Sexual Maturation 

Control of puberty is an important issue in fish farming. 
Precocious sexual maturation has economic implications, 
and various methods have been devised to attain the efficient 
production and public acceptance in terms of welfare and 
sustainability (Taranger et al. 2010). Selective breeding, envi- 
ronmental manipulation, induced triploidy, and monosex pro- 
duction are among the methods used to control the sexual 
maturation. The brain-pituitary-gonad axis regulates sexual 
maturation in teleosts. Precise mechanisms are not fully 
understood. Several factors influence sexual maturation, in- 
cluding habitat, temperature, photoperiod, nutritional status, 
social interaction, pheromones, and hormones, all having role 
in the extensive regulatory feedback system (Taranger et al. 
2010). 

Little is known about the role of miRNAs during sexual 
maturation in teleosts. miRNA profiling studies in Atlantic hal- 
ibut showed significant differences in the expression of many 
miRNAs. The expression of let-7a, miR-1 43, miR-1 45, and 
miR-202-3p was significantly higher in adult testis compared 
with adult ovary, and miR-451 was significantly downregu- 
lated in brain of juveniles compared with adult females 



1926 Genome Biol. Evol. 6(8): 191 1-1937. doi:10.1093/gbe/evu151 Advance Access publication July 22, 2014 



miRNA in Teleost Fish 



GBE 



(Bizuayehu et al. 201 2b). In Nile tilapia, the expression of miR- 
129-3p and miR-727-3p was significantly higher in mature 
females than males, whereas the expression of miR-132a 
and miR-212 was significantly higher in mature males than 
females (Xiao et al. 2014). Differential expression of miRNAs 
between mature and immature gonads has been reported in 
rat (Gaytan et al. 2013), pig (Luo et al. 2010), and chicken 
(Kang et al. 2013). Change in expression of miR-1 45/c-Myc/ 
Lin-28/let-7 axis in hypothalamus occurs during rat puberty 
(Gaytan et al. 2013). In rainbow trout, 13 miRNAs showed 
differential expression patterns during the ovarian develop- 
ment (Juanchich et al. 201 3). Further research on mechanisms 
involved in the regulation of sexual maturation through epi- 
genetics, and miRNA-related pathways would advance our 
understanding of fish sexual development. 

miRNA in Endocrine Organs of Teleosts 

Hormone-secreting organs in teleosts include the brain 
(hypothalamus, pituitary, and pineal gland), thyroid, kidney 
(chromaffin tissue and corpuscles of stannous), gonad 
(theca and Leydig cells), intestinal mucosa, pancreatic islets, 
ultimobranchial body, and urophysis. The development and 
tissue specification of these organs involve miRNAs. 
Expression profiling studies in fish indicate tissue specificity 
of miRNAs in some of these organs (Wienholds et al. 2005; 
Kapsimali et al. 2007; Kloosterman et al. 2007; Tessmar-Raible 
et al. 2007; Bizuayehu et al. 2012b). 

Brain 

Endocrine compartments of the brain: Hypothalamus, pitui- 
tary, and pineal gland are involved in the regulation of 
homeostasis. 

The Hypothalamus 

The hypothalamus comprised preopticus, lateralis tuberis, 
recessus lateralis, and recessus posterioris nuclei, which pene- 
trate the pituitary to promote hormonal releases (Goos 1 978). 
Few studies show expression of miRNAs in hypothalamus of 
teleosts (Wienholds et al. 2005; Tessmar-Raible et al. 2007). 
Some of these miRNAs, such as miR-7 and miR-7b, are found 
both in mice (Bak et al. 2008) and zebrafish (Tessmar-Raible 
et al. 2007), indicating probable functional conservation. 
Although there is no report available on miRNA expression 
in different nuclei in teleost hypothalamus, miRNA profiling 
of mammalian homologous regions indicates differential 
and specific expression patterns (Herzer et al. 2012). 
Corticotrophin-releasing factor from hypothalamus enhances 
proopiomelanocortin transcription and promotes adrenocorti- 
cotropic hormone production. This process is negatively reg- 
ulated by miR-375 in mice (Zhang et al. 2013a). Another 
corticotrophin-releasing factor, urocortin 2, promotes expres- 
sion of miR-325-3p in pituitary and suppresses biosynthesis 



and secretion of luteinizing hormone in rat (Nemoto et al. 

2012) . Gonadotropin-releasing hormone regulates multiple 
miRNA expression in gonadotrope cell lines, producing 
downregulation of miR-99b and miR-125b, and upregulation 
of miR-132, miR-151, miR-212, miR-222, miR-350, and miR- 
424 (Godoy et al. 201 1 ). Intronic miR-1 32 and miR-21 2 target 
p250RhoGAP, thereby guide morphological change and 
increase the motility of gonadotropes (Godoy et al. 2011). 
In chicken, hypothalamic neuronal cell migration is mediated 
through miR-1 38 repression of RELN (Kisliouk and Meiri 

2013) . 

The Pituitary Gland 

One of the functions of the pituitary gland is signal transmis- 
sion between hypothalamus and peripheral tissues. In tele- 
osts, pituitary is composed of neurohypophysis, a neural 
component originated from the diencephalon, and adenohy- 
pophysis derived from the buccal epithelium (Schreibman 
et al. 1973). Adenohypophysis is divided into rostral pars 
distalis or proadenohypophysis, proximal pars distalis or 
meso-adenohypophysis, and pars intermedia or meta- 
adenohypophysis (Schreibman et al. 1973). There are many 
signaling molecules and transcriptional factors that control 
teleost pituitary development and patterning (Sbrogna et al. 
2003; Herzog et al. 2004; Nica et al. 2006). Some of these 
transcriptional factors are known to be targeted by miRNAs in 
mammals; for example, miR-26b represses lymphoid en- 
hancer factor 1 to promote the generation of somatotrope, 
lactotrope, and thyrotrope cell lineages in mice (Zhang et al. 
2010). In contrast to mammals, there are interdigitations 
between neurohypophysis and adenohypophysis in teleosts. 
Considering the role of miRNAs in defining boundaries, such 
as midhind brain boundary (Leucht et al. 2008) and hindbrain 
and spinal cord boundary (Woltering and Durston 2008), the 
presence of functional distinction in amalgamated tissues 
raises the question whether miRNAs have a role in demarca- 
tion of the pituitary morphological boundaries. Teleosts rostral 
pars distalis is composed of various cell types, such as r|, e, 
neck, and channel cells (Schreibman et al. 1 973). r| and e cells 
produce prolactin and adrenocorticotrophic hormone, respec- 
tively, whereas other cells are involved in the movement of 
material into or from the follicular lumina, as well as in phago- 
cytosis and dispersal of a hormone or a carrier substance into 
peripheral circulation (Schreibman et al. 1973; Leatherland 
1976). Teleost proximal pars distalis is composed of cells se- 
creting thyrotropins, somatotropins and gonadotropins, 
whereas pars intermedia contains two cell types, which se- 
crete melanocyte stimulating hormone and somatolactin 
(Rand-Weaver et al. 1991; Schreibman et al. 1973). 
Differentiation and functional specification of these cells re- 
quire specific regulatory elements promoting (Liu et al. 2003) 
or inhibiting (Liu et al. 2006) lineage-specific gene expression. 
The role of miRNAs in such the specification and in release of 
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neurosecretory material remains unknown. In mammals, miR- 
7, miR-7b, miR-141, miR-200a, and miR-375 are enriched in 
the pituitary (Landgraf et al. 2007; Bak et al. 2008); similarly, 
in situ detection during zebrafish development shows the 
expression of miR-375 in pituitary (Wienholds et al. 2005; 
Kapsimali et al. 2007), indicating evolutionary conservation 
of miRNA in pituitary function. 

The Pineal Gland 

The pineal gland modulates physiological activities related to 
daily and seasonal rhythmicity in fish. Production of the main 
compound, melatonin, is stimulated by darkness and inhibited 
by light. Involvement of miRNAs in teleost pineal gland devel- 
opment and secretion is unknown. In rat pineal gland, expres- 
sion of several dominant miRNAs such as miR-96, miR-124, 
miR-125b, miR-127, miR-182, and miR-183 has been found 
(Clokie et al. 201 2). Although majority of miRNAs identified by 
the authors had similar expression during the day and night, 
the 2-fold increase in miR-96, miR-182, and miR-183 during 
the day when compared with the night was found. Moreover, 
a significant reduction of melatonin in miR-483-transfected 
pinealocytes and repression of arylalkylamine A/-acetyltransfer- 
ase expression, an enzyme converting serotonin to N-acetyl- 
serotonin, was observed. It is essential to consider the 
differences in photoperiodic and circadian control of neuro- 
endocrine functions between teleosts and mammals. The 
mammalian system has a linear flow (eye-hypothalamus- 
pineal gland) to the rhythmic production of melatonin, 
whereas in fish, melatonin biosynthesis and degradation re- 
quire synchronized complex interactions among eye, pineal 
gland, brain (particularly hypothalamus and pituitary), and pe- 
ripheral tissues upon the stimulation with light (Falcon et al. 
2010). Therefore, functional characterization of miRNAs in 
teleost pineal gland will be important to understand neuroen- 
docrine regulation of rhythmicity by the pineal gland. 

The Thyroid Gland 

The thyroid gland secretions, thyroid hormones (THs), are con- 
served among vertebrates. They mediate gene expression by 
binding to thyroid-hormone-binding protein, which binds to 
TH responsive element and regulates the transcription of dif- 
ferent genes (Wu and Koenig 2000). In teleosts, THs are ma- 
ternally deposited in eggs (Brooks et al. 1997). They are 
essential for growth and metabolism throughout the life 
(Power et al. 2001 ; Liu and Chan 2002) and regulate different 
aspects of development including metamorphosis, which 
requires substantial changes in morphology, physiology, and 
behavior (Yamano et al. 1991; Liu and Chan 2002). 
Metamorphosis is accompanied by alteration in miRNA profile 
(Fu et al. 201 1, 2013; Bizuayehu et al. 2012a). In rat thyroid 
cells, many miRNAs such as miR-1, miR-28a, and miR-296-3p 
are differentially expressed and possibly they target transcripts 
that are important in thyroid cell proliferation (Leone et al. 



201 1). Akama et al. (2012) demonstrated that the expression 
of 47 miRNAs in rat thyroid cells was reduced after the addi- 
tion of thyroid stimulating hormone. Their result suggests that 
thyroid stimulating hormone regulates thyroid cell prolifera- 
tion partly by reducing the expression of miR-1 6 and miR-1 95, 
which target important cell proliferation genes, including 
Mapk8, Ccnel, and Cdc6. Similarly, using Dicer conditional 
knock-out model, Frezzetti et al. (2011) showed the impor- 
tance of miRNAs in the morphology and function of thyroid 
gland in mice. 

The Endocrine Pancreas 

The endocrine pancreas physiology has been studied in fishes 
for more than a century. Pancreatic hormones: Insulin, gluca- 
gon, glucagon-like peptide, somatostatin, pancreastatin, and 
pancreatic peptide have been isolated and localized in differ- 
ent cells of the pancreas (Plisetskaya 1 989; Jonsson 1 991 ). The 
secretion of pancreatic hormones is a tightly controlled pro- 
cess with feedback mechanisms. Little is known about the 
physiological function of miRNAs in teleost pancreas 
(Wienholds et al. 2005; Kloosterman et al. 2007). Some 
miRNAs, such as miR-375, are uniquely expressed in pancre- 
atic cells; Kloosterman et al. (2007) have shown the impor- 
tance of miR-375 in the insulin-secreting pancreatic islets, as 
miR-375-knockdown zebrafish embryos had dispersed islet 
cells. In mammals, miR-7, miR-9, miR-29b, miR-30d, miR- 
124a, and miR-375 regulate the secretion and islet develop- 
ment (Poy et al. 2004; Baroukh and Van Obberghen 2009; 
Tang et al. 2009; Pullen et al. 201 1). Poy et al. (2004) have 
demonstrated that miR-375 targets myotrophin and is in- 
volved in insulin exocytosis. Similarly, Plaisance et al. (2006) 
have reported the control of the secretory function of insulin- 
producing cells by miR-9. In their electrophoretic mobility shift 
assay, chromatine immunoprecipitation, and gene reporter 
experiments, the transcriptional factor onecut-2, which tar- 
gets granuphilin, is implicated in insulin secretion, and its 
level is kept at an appropriate level by miR-9. Given a persis- 
tent hyperglycemia in some aquaculture species fed carbohy- 
drate-rich diets (Moon 2001), studies on carbohydrate 
metabolism in the context of gene regulation may help in 
understanding the consequences of fish meal replacement 
with plant products in aquaculture feeds. Therefore, a catalog 
of abundant and specific miRNAs in pancreatic cells of teleosts 
and comparison of expression patterns and types of miRNAs 
among different phylogenetic groups will advance our under- 
standing of endocrine pancreas in teleosts. 

Other Endocrine Cells/Tissues 

There are no reports yet on miRNA characterization in other 
endocrine tissues of teleosts, including chromaffin tissue, cor- 
tical tissue, hormone secretory cells in gonad, intestinal 
mucosa, or urophysis. 
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Response to Environment 

Osmoregulation 

Osmoregulation in teleosts is a process whereby an organism 
adapts to different ionic environments. In zebrafish embryos, 
miR-8 family (miR-8, miR-141, miR-200a, miR-200b, miR- 
200c, and miR-429) has been expressed in ionocytes (Flynt 
et al. 2009), which are specialized branchial epithelial cells 
that are involved in the maintenance of osmotic homeostasis. 
An increase in osmotic stress sensitivity of miR-200a and miR- 
200b morphants was observed along with Na + accumulation 
in ionocytes, indicating the function of miR-8 family in Na7H + 
exchanger (2009). In adult Nile tilapia, apart from high expres- 
sion of miR-30 family and miR-429 in kidney and gills, miR-30c 
and miR-429 are implicated in osmotic stress regulation (Yan, 
Zhao, et al. 2012; Yan et al. 2012b). It is not known whether 
the different cells that are involved in osmoregulation respond 
in similar fashion, what is celLmiRNAs specificity, and whether 
similar miRNAs are involved in osmoregulation in marine and 
freshwater species. Participation of miRNAs in osmotic stress 
regulation is disclosed through the modulation of genes that 
are involved in regulation of membrane dynamics, trafficking 
of transmembrane proteins, and osmoregulatory signaling, 
such as Na + /H + exchanger regulatory factor 1 in zebrafish 
(Flynt et al. 2009) and osmotic stress transcriptional factor 1 
in Nile tilapia (Yan, Zhao, etal. 2012). Considering functions of 
different hormones, such as growth hormone, somatolactin 
and prolactin in ion control in teleost fish (Rand-Weaver et al. 
1991), a large regulatory pathway model can be anticipated. 
Elucidating miRNA functions in passive and active transporters 
through the regulation of membrane-bound proteins, periph- 
eral proteins, and/or porin structures would help to under- 
stand osmoregulatory mechanisms. 

Temperature 

Temperature is a major external determinant of teleost phys- 
iology. All teleosts have a preferred water temperature, and a 
deviation from the optimum range can affect survival, growth, 
behavior, and reproduction. miRNAs have important role in 
regulation of environmental effects (Li et al. 2009); however, 
little is known about their role in teleosts thermal regulation. 
Although 11 upregulated and 15 downregulated miRNAs 
were found in brain of cold-acclimated zebrafish, their overall 
contribution in the regulation of protein-coding genes was 
found to be minimal (Yang et al. 2011). The authors sug- 
gested that brain miRNAs had developmental rather than 
thermal adaptation functions. More work is needed to test 
this hypothesis using different species, tissues, and experimen- 
tal set-ups. Significant differences in miRNA expression be- 
tween hyperplasic myotube and hypertrophic phenotypes 
have been found in zebrafish as the effect of temperature in 
early development (Johnston et al. 2009). This report also 
showed significant differences in myotomal fast muscle fiber 



recruitment in embryos held at different incubation tempera- 
tures (22 °C, 26 °C, and 31 °C). Thus, it is interesting to know 
the role of miRNA in this process. Further work on both eu- 
rythermal and stenothermal teleosts is needed to understand 
the functions of miRNAs in thermal regulation. It is important 
in the context of global warming in general and rise in oceans 
temperature in particular. Knowledge on the miRNA regula- 
tory modulation as a function of temperature can be useful to 
model ecological and biological consequences of the global 
warming. 

Stress 

Transcriptional and posttranscriptional adjustments of gene 
expression during the stress are among the long-term re- 
sponses, and they are primarily regulated by hormones (Fiol 
and Kultz 2007). The accumulation of Argonaute proteins and 
miRNAs in separate cytoplasmic foci in stress-dependent 
manner (Leung et al. 2006) indicates possible modulation of 
stress response by miRNAs. Different stressors, such as oxida- 
tive, osmotic, or heat shock can alter miRNA expression. 

Oxidative Stress 

Oxidative stress is caused by reactive oxygen species and has a 
potential to inhibit cellular functions, but natural antioxidants 
such as tocopherols (vitamin E) can limit the damage. Feeding 
Nile tilapia a diet containing various levels of vitamin E (0, 50, 
and 2,500 mg/kg) resulted in reduction in the expression of 
miR-16, miR-122, miR-146a and miR-223 in E-deficient group 
and increased the expression of miR-16, miR-21, miR-122, 
miR-125b, miR-146a, miR-155, miR-181a, and miR-223 in 
E-enriched group, when compared with E-normal group 
(Tang et al. 2013). Yu et al. (2010) have shown antioxidant 
activity of miR-451 in zebrafish embryos. Further investiga- 
tions on the role of miRNAs in oxidative stress and their po- 
tential as biomarkers are important for aquaculture to improve 
welfare and growth performance of fish. 

Osmotic Stress 

Uptake of hyperosmotic 2% saline water resulted in upregu- 
lation of expression of miR-7b, miR-9, miR-29b, miR-1 37, and 
miR-451 and downregulation of miR-409, miR-1 07, miR-1 03, 
miR-1 85, and miR-320 in hypothalamus in mice (Lee et al. 
2006). This study showed also a reduction of Fos protein pro- 
duction (component of a transcriptional factor activator pro- 
tein 1) in paraventricular nucleus as the result of miR-7b 
activity. Similarly, salt tolerance of Nile tilapia is partly attrib- 
uted to regulatory function of some miRNAs, such as miR-30c 
and miR-429 (Yan, Zhao, et al. 2012; Yan et al. 2012b). 
Nevertheless, the mechanism of miRNA action in osmotic 
stress response is not known yet. 
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Irradiation, Chemical, and Physical Stressors 

In zebrafish embryos, miR-125b is downregulated by gamma- 
irradiation and upon exposure to a cytotoxic chemical, which 
corresponds to an increase in p53 protein level; it indicates the 
possible function of miR-125b in stress response regulation 
(Le et al. 2009). Physical stress can act as a "morphogen"; 
for example, shear stress and stretch of cardiomyocytes acti- 
vate the differential expression of miR-21 during heart valve 
formation in zebrafish (Banjo et al. 2013). Exposure of lake 
whitefish (Coregonus lavaretus) to microcystin significantly 
changed the expression of 6 liver miRNAs in a time-depen- 
dent manner (Brzuzan et al. 201 2), indicating the involvement 
of miRNAs at different levels of physiological acclimation re- 
sponses. Likewise, 24 h exposure of adult zebrafish to differ- 
ent concentrations of microcystins altered the expression of 
4 miRNAs (Li et al. 2013). Although the extract, delivery 
method, and concentration differed in those two reports, 
miR-122 expression pattern was opposite in whitefish and 
zebrafish liver at the similar exposure time. Exposure 
to other chemicals, such as T2,3,7,8-tetrachlorodibenzo-p- 
dioxin and perfluorooctane sulfonate also altered the expres- 
sion of group of miRNAs during the zebrafish development 
(Zhang et al. 201 1; Jenny et al. 2012). All these studies indi- 
cate that miRNA expression is affected by a variety of environ- 
mental stressors; however, the mechanisms of this 
modulation and their consequences need to be further 
elucidated. 

In Vivo Models for miRNA Studies 
in Teleosts 

In vivo models are useful tools for expanding knowledge 
about the molecular mechanisms. They are not widely used 
yet in miRNA research in teleosts. 

The first teleosts cell line, RTG2, has been derived from a 
gonad of rainbow trout (Wolf and Quimby 1962). Since then a 
number of fish cell lines have been established and utilized 
in many fields of biological research. A rainbow trout spleen 
cell line, RTS34, was used in target validation of a miRNA in 
Atlantic halibut, which confirmed the binding of miR-24 to the 
3'-UTR of kiss peptin 1 receptor-2 (Bizuayehu et al. 2013). 

Apart from in vitro studies, miRNA mutant and transgenic 
lines should be explored as these lines in zebrafish have wid- 
ened our understanding of the molecular processes that gov- 
ern the phenotypic outputs of several transcriptional factors, 
signaling molecules, and genes. However, functional versatility 
of miRNAs during the transcription and posttranscription pro- 
cesses suggest careful reexamination of the existing knock- 
out/mutants models in terms of miRNA function. 

Recent genome editing technologies, such as zinc finger 
nucleases, transcription activator-like effector nucleases 
(TALENs), CRISP-Cas, or RNA-guided Cas9 nuclease (Urnov 
et al. 2010; Bedell et al. 2012; Chang et al. 2013; Cong 



et al. 2013; Hwang et al. 2013) provide an opportunity to 
edit a specific miRNA or a cluster in the genome. For example, 
two miRNA clusters, miR-17-92 and miR-430, were success- 
fully interrupted in zebrafish using TALEN technology (Liu et al. 
2013). The functional characterization of miRNAs will be ac- 
celerated by integrating these methods with spatiotemporal 
miRNA transcript regulation (knock-down and overexpression) 
using various methods such as caged antagomirs (Connelly 
et al. 2012; Griepenburg et al. 2013). The prospect in this 
regard is immense. 

Conclusions 

Functions of miRNAs in teleost development, growth, and 
physiology are not satisfactory understood yet. Evolutionary 
diversity of teleosts and the resulting plasticity in their environ- 
mental adaptation suggest that there is much more to dis- 
cover beyond the very few species and processes investigated 
so far. This review indicates the importance of functional stud- 
ies. Although a given miRNA may be either predicted (in silico), 
or demonstrated (in situ), to have a regulatory capacity, its 
physiological relevance has to be established in a specific con- 
text of a live system. Characterization and profiling of miRNAs 
in different teleosts is the beginning of a long road ahead to 
disclose the functional conservation and variations of miRNA 
regulatory pathways across Teleostei. 
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